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We show the possible stable soliton generation for the dark-state polaritons (DSPs) in an 
electromagnetic induced transparency (EIT) medium composed of A-type atoms. Whether the 
solitons are dark or bright can be controlled by the coupling field intensity and the one photon 
detuning of the probe field. The velocity, spatial and time widths of the solitons can also be 
adjusted by the coupling light. 
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The intriguing problem to find new quantum systems whose wave functions can form solitons and 
to study their novel dynamical properties always attracts considerable interests. The well-known 
examples include the Ginzburg-Pitaevskii-Gross equation in Bose-Einstein Condensate(BEC) 

and the Maxwell-Bloch equation in nonlinear electric medium, where there exists soliton solutions 
that respectively describe the properties of wave functions of the atomic condensate and the electric 
field. A natural question then may be asked: in an interacting system of atoms and electromagnetic 
field, if we treat them as a total system, can the coupled matter-photon system also form a soliton? 
An important concrete example may be the light storage in the electromagnetic induced transparency 



(EIT)P] medium which, in recent years, received much attention due to its potential applications in 
the field of quantum information science. 

In fact since the technique of resonant enhancement of the index of refraction without absorption 
was proposed P] and many accompanying striking phenomena was observed [EflSSEjQl i the light 
storage with the technique of EIT has been an exciting research field in current literature, especi ally 
after the "dark-state polaritons" (DSPs) theory was proposed by Fleischhauer and Lukin et al. jl2Lll,l | 
based on a field theory reformulation of the adiabatic approximation |^^. DSP is a new quantum field 
which is the superpositions of the electric field amplitude e{z, t) and atom coherence pcb between two 
lower levels of the A type atoms, and it describes the total system of the electric field and collective 
atomic state. In linear theory where the one-photon detuning of the probe pulse is zero, the quantum 
state of the polaritons can be mapped from the electric field state into the collective atomic excitation 
without change when the coupling laser was adiabatically turned off Q| , which implies that DSP 
is really an elegant description of the total system. 

The paper will discuss the nonlinear properties of the DSPs and prove that the motion of the DSPs 
satisfies a (l+l)dimcnsional nonlinear Schrodinger equation (NLSE) that has possible stable soliton 
solutions in an electromagnetic induced transparency (EIT) medium composed of A type atoms. Note 
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that the soliton is not formed by the probe pulse or atom coherence alone, but by the total state 
function of them. This is different from the familiar conventional ones|lfj| which are formed by the 
wave function of a single physical system. The coupling field intensity, along with the one photon 
detuning of the probe field, is shown to decide whether the solitons are dark or bright and other 
important parameters, like the time width and the velocity of the solitons. 

We consider the quasi 1-dimensional system composed of three-level A-type configuration atoms 
with energy levels assumed to he Ea > Ec > Ei,. A coherent probe field with positive frequency part 
of the electric field Ep'^\z,t) couples the transition between the ground state |&> and the excited 
state \a>. A = ujab — lo is the small one-photon detuning between the carrier frequency u and the 
atomic transition frequency ujat- The stable state |c> is coupled to \a> via a coherent coupling field 
with Rabi- frequency Q. We assume the coupling pulse is much stronger than the probe one and its 
frequency is ly = ujac- Due to the small intensity of the probe field, we expand the density matrix in 
the following form: 

where /i, = a, b. The quantities p^^} is of the same order of smallness as intensity of the probe 

(2) 

pulse, the pin) is of the second order of the smallness, and so on. To analyze the nonlinearity of 
the susceptibility of the probe light, we calculate the density matrix element pat to the third order. 
Together with the relation: pafc"a(pift^ ~^pfb) — ^o(x'^^ {^)+X^^^ where £ is the dimensionless 

slowly varying amplitude of the probe pulse Ep{z, t) and Ua being the atom density, the susceptibility 
(x ~ x'^^ + 3x^'^''kP) is then given by 

where ^^i, {p, v = a,b, c) is the transverse decay rate between levels p and v and here they satisfy: 
7ch,7ca 9 — Pab^J^^f^^ , and N is the effective total atoms in the quantum volume V. The 

susceptibility is related to the refractive index through n = -^1 + XC^); then the wave vector of the 
probe pulse can be approximately calculated by 

3 

kKkab+ 2fcabX^^^kP +(il{uJ~ UJab) + /32(t^ - ^^abf (3) 

where /3i = l/(ccos^ is the inverse of the group velocity of the probe pulse, (32 = 

5f2^|w=w„6 = — tan'' 0/(a;oc) and fcaj, = oJab/c. The mixing angle 6* is defined via tan 6'(t) = g\fN /VL{t), 
which is a sufficiently slowly time-dependent function. Assuming that the wave vector of the probe 
pulse has a narrow spreading around the central value fco: from the formula (3) the frequency of the 
probe pulse can be expanded around fco in the following form: 



^0 - TT^oK^X^'^kl' + Vg{k - ko) - P2Vhk - kof (4) 



where Vg — ccos^ 9 and ujq ~ k^c are respectively the group velocity and central frequency of the 
probe pulse. So the dispersion equation of probe field Ep{z,t) can be approximately given by 

^dEp^^ (c.o-ifcoK,x(')kni?p(^,<) +Vg{~^-ko)Epiz,t) 

~P2VH\f^-kofEp{z,t) (5) 
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In above derivation the varying of Rabi frequency with time is ignored. In fact, the time- 
dependent character of brings a reduction/enhancement of the amphtude e{z,t) by a contribution 
2(i+a) ^ ^^^^to the motion equation of the probe pulse during the Raman adiabatic passage, where 
the group index a — g^N/OF' and a — da/dt. Consequently the propagation equation of the probe 
pulse in the EIT medium is 

\ d d OL 

*fco(-^^ + ^ + ^)£(2,i) = k^hV^^e{z,t) -^{LOQ)\e\^£{z,t) (6) 

(3) 

where ?7(wo) = — fc^ — ^'^^ vacuum light speed. 

As is known that the propagation of the electromagnetic pulse in an EIT medium can be easily 
understood in terms of new quantum fields, i.e., polariton-like fields, which are superpositions of the 
dimensionless electric field amplitude s(z^t) and the atom coherence p^h- The dark-state polaritons 
and bright-state polaritons(BSP) here are defined byQ 



^{z,t) = cos9{t)e{z,t) - sin9{t)V N pcb{z,t) exp{i{kQ - kc)z) (7) 



$(2;,t) = sm9{t)e{z,t) + cos 9 {t)vNpcb{z,t) ex.p{i{kQ - kc)z) (8) 

where kc is the wave- vector of the coupling field in z direction, the mixing angle ^ is a function of 
time, while the electric field e(z, t) and atom coherence pcb are functions of time and z-axis coordinate. 
Together with the formula@ , by a straightforward calculation one can derive the propagation equation 
of the DSPs: 

'koi^m + - fco/32if ^*(^, t) 

= ^ko^(>^ - ^ sin0(^ + Vn^) ~ zfcotan0^$ 

+kop2VgHan9-§^^ - 7^iiOo)\£\''e{z,t) (9) 

where Pcb = Pcb{z,t) exp(i(fco — kc)z) and 9 = d9/dt. The motion equation of BSPs satisfies [l^ : 

^ = Pw{wt+^ba + i^) (tan 0) [(f + zA)(sin * - cos $) + g^epca] + * f^^ba- (10) 

where Fba is Langevin noise term which will be omitted in the following derivation; and the density 
matrix \pca\ ^ and \pca\ ~ \^^Pcb\ \pcb\ due to the weak probe field and low atomic excitation. 
By introducing the adiabatic parameter e = {g^/NT)~^ with T being a characteristic time 
we find in the lowest order adiabatic theory that, |<i>(z,<)| ^ IffPcal ~ jf^lff Pi cos6'^'| for the 
ultraslow light case. It is clear that when jfip > |f;A| or jfip ^ Is^L one has |<i>| ^ |cos0^| 
(since the control field is much stronger than the probe one, i.e., jfip ^ Iff^P)- The typical values 
of these parameters in ultraslow light case are [l^ [l^ Q. = 1.0 x 10^~*s~^, A — 1.0 x 10^s~^ 
g^/N = 5.0 X IQi^s-i and 5 = 2.0 x 10^ from which we find cos9 - lO^^-^-^) and - \gA\ and 
then |<I>| ^ lUpjcosS^I ^ |cos0^|. Hence one can approximately obtain £{z,t) k, cos 9 {t)'i/{z,t), 
Npcb ~ — sin 9 {t)'^(z,t) and $ w 0. Substituting these results into the formula JHl yields 

ik^,^{^,t) + g,^{tt) - -C,,\^\'^{tt) (11) 
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where k = k^) / {Vg sm^ 9) , the slowly varying nonlinear coefficient C„ = ri{uj^)) co\? csc^ 9 , and ^ = 
z — J^^ drVg is the coordinate in the rest frame of the probe pulse. This is a (l+l)-diniensional 
nonlinear Schrodinger equation (NLSE), which has a bright (dark) soliton solution when Ke{ri{u!o)) > 
{Re{i]{uJo)) < 0). As is known, the NLSE with slowly varying coefficients can be solved with 
perturbation theory 18]. The substitution t' = J^dt/k, ^'(^,i) = transforms Ea. (|ll|l into 



standard form with perturbation *^2^^^^ llaUi 

lir^'i^, t') + ^vl;'(^, t') + 2\^r^'{t t') = lePm (12) 

whereP(^'') = 2C~^I^*'(^' ultraslow light case Vg < c, we have ^(vE') w ianO"^ . Then, 
when A < and <^ < il^, with a perturbation theory and one can obtain the fundamental 

mm 



bright soliton 

rt' 



exp(i / dT2A ) 

Jto 

where A = cos 0(0)/ cos 9{t) and Aq is a constant related to the initial condition. Then 



(13) 



*b = A/sec/i(VlCVjAi"Oexp(W dTM'^\Cn\/k) (14) 



or 



to 



t rt 

2| 



^-fc = Afsec/i(v1ajM(z - / VgdT))exp{i dTM^\Cn\/k) (15) 

where M = (t)/2|°/2 rose(t) maximum amplitude of the soliton. From the formula 

(|13|l and (|14|l one can easily find the spatial width of the soliton A^ = ln{2 + ^/3)\Nc{Q^ — 
A'^)/{6kon'^AP{t)A)\^/'^ cos9{t) oc fl{t), which can be easily controlled by the coupling light. When 
the Rabi frequency is adiabatically reduced, the spatial width narrows according to A^(t) = 
A^(0)r2(t)/ri(0). Likewise we obtain the time width of the soliton tpwhm ~ 

ln{2 + V3)\N{n^ - 

A^)/{6kofPM^{t)A)cos'^ 9{t)\^/^ cx l/fl{t), which is inversely proportional to the Rabi-frequency 
of the coupling light. When is adiabatically reduced, the time width broadens according to 
TpwHAiit) — TpwHMiiO)^iO)/^it). The result that time and spatial widths can be controlled by the 
coupling light is an important observation for the properties of the bright solitions formed by DSPs. 
To give a numerical estimate of the spatial width of the soliton, we assume the Rabi frequency 

of the coupling pulse = 1.0 x 10®s"\one-photon detuning A = 1.0 x 10'^s"\ g^/N = 5.0 x 10^^s"\ 

= 1.0 X 10^^, AqCOs9{0) w 0.1, and ujq = 2.0 x 10^^s~\ the spatial width can be then calculated 
A^ « 2.63 X lO^^mm. This estimation indicates that this interesting result may be observable under 
present experimental techniques (the length of the EIT medium used in present experiment is about 
several centimeters [21] ) . By choosing Eq — M(0), L = I/Aq, and r — g'^N/{Q'^{Q)Aoc) as units to 
normalize the amplitude of the DSPs, coordinate z and time t, respectively, the amplitude of the 
fundamental bright soliton can be plotted in Fig.l. 

The derivation of the above result is under the condition Re(?7(cjo)) > 0. However, when the 
coupling field intensity is adiabatically reduced to fl^ < A^ while even il^ ~ |(;A| with A < 0, the 
nonlinear coefficient C„ becomes a negative. Meanwhile the NLSE supports a dark soliton solution. 
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the amplitude of which is a hyperbohc tangent function. The time or spatial width of the dark solitons 
has the same form as that of the bright ones. On the other hand, if we initially set a positive one- 
photon detuning, i.e. A > 0, the DSPs can also form dark solitons (for the case fi^ > A^) and bright 
ones (for the case < A^). 

It is noteworthy that in our formulation we have ignored the transverse decay rate jab which can 
lead to a contribution — i(5|^'p^(z, i) to the right-hand side (RHS) of eq.(12), where the coefficient 

1 /2 

S cx Im(C„)/|C„| w 7afc|A|/[(ri^ — A^)^ + 7^f,A^] represents a nonlinear loss (absorption) and 
instability for the propagation of NS solitons. For present purpose, when — A^| 3> |7ahA|, the 
effect of nonlinear absorption can be ignored. The velocity of the solitons we discussed above is just 
the group velocity of the probe pulse and it can also be conveniently controlled by the coupling pulse. 
Another interesting issue is that we may extrapolate the idea in this paper to other interacting systems 
such as a passive medium j22| or a Tonks gas [2,'^ . and we may even introduce a new technique for 
quantum memory with solitons formed by total state function of quantum memories (QMEs) and 
quantum carriers (QCAs). Due to the strong stability of solitons, quantum information may be 
elegantly stored and released against some environmental noise or external field disturbances in this 
technique. 

In conclusion we have shown that the dark-state polaritons in an EIT medium can form the possible 
stable solitons, and whether the solitons are dark or bright is dependent on the coupling field intensity 
and the one-photon detuning of the probe field. Note that the kind of solitons arc formed by the total 
wave function of two different physical systems, i.e., the electromagnetic field and collective atomic 
excitations (spin waves). The velocity, time and spatial widths of the solitons are also shown to be 
adjustable by the coupling light. Of course, these results are derived for a weak probe field and low 
atomic excitation, i.e., the atomic density is approximately be treated as constant, and the more 
interesting case with a varying density is not studied at present. This challenging issue, together with 
the evolution process from bright solitons to dark ones, deserves further studies in our future works. 
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FIG. 1: The amplitude of the fundamental bright soliton formed by the DSPs with dimensionless variables. 
The normalized factors are Eo = M(0), L = 1/Ao, and t = g'^ N/(Jl^ (0)Aqc) that represent the units of the 
amplitude of the DSPs, coordinate z and time t. 



